function y = computeDamageFunction(element,etilde) %compute F(etilde)
Dnn = element.Dnn;
g = computeTractionLaw(element,etilde);
y = Dnn*etilde/g;
end
